Time Delay Induced Death in Coupled Limit Cycle Oscillators 
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We investigate the dynamical behaviour of two limit cycle oscillators that interact with each other 
via time delayed coupling and find that time delay can lead to amplitude death of the oscillators 
even if they have the same frequency. We demonstrate that this novel regime of amplitude "death" 
also exists for large collections of coupled identical oscillators and provide quantitative measures of 
this death region in the parameter space of coupling strength and time delay. Its implication for 
certain biological and physical applications is also pointed out. 
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Coupled limit cycle oscillators provide a simple but 
powerful mathematical model for simulating the collec- 
tive behaviour of a wide variety of systems that are of 
interest in physics |p]-p^, chemistry and biologi- 

cal sciences . These oscillators have also attracted 

some large scale numerical [ p^ and novel experimental ef- 
forts [|l6| . For weakly coupled oscillators the predominant 
effect is a synchronization of the frequencies of the indi- 
vidual oscillators to a single common frequency once the 
coupling strength exceeds a certain threshold, while the 
amplitudes remain unaffected. For stronger couplings the 
amplitudes also play an important role and give rise to 
interesting phenomena like the Bar-Eli effect ||l^ where 
all the oscillators suffer an amplitude quenching or death 
. In general there can be a wide variety of collective 
behaviour including partial synchronisation, phase trap- 
ping, large amplitude Hopf oscillations and even chaotic 
behaviour |^|^. In recent times there have been exten- 
sive investigations of coupled oscillator systems including 
elegant statistical mechanics formulations in the limit of 
infinite number of oscillators M.K(M . 



The salient features of the behaviour of a finitely large 
number of oscillators (usually obtained from numerical or 
approximate analytic means) can often be understood by 
analysing just two coupled oscillators. We have carried 
out such an analysis to investigate the effect of time de- 
lay on the interaction between two limit cycle oscillators. 
Time delay is ubiquitous in most physical and biological 
systems p^-n9[, arising from finite propagation speeds 
of signals for example, and have not been widely studied 
in the context of coupled limit cycle oscillator systems. 
Niebuhr et al and Schuster and Wagner who are 
one of the few who have carried out such an investigation, 
have restricted themselves to the simpler coupled phase 
models where the phenomenon of amplitude death does 
not exist. In our model equations we have retained both 
the phase and amplitude response of the oscillators and 
we find that time delay has a significant effect on the 
characteristics of all the major cooperative phenomena 
like frequency locking, phase drift and amplitude deaths. 
In particular our detailed numerical investigations show 
that in the presence of time delay the parameter regime 



of amplitude death can extend down to the region of 
zero frequency mismatch between the oscillators. This 
is in sharp contrast to the situation with no time de- 
lay where all previous numerical and analytical studies 
lH'lljlJll show that amplitude death can occur only if 
the coupling between oscillators is sufficiently strong and 
when the frequencies are sufficiently disparate. In this 
Letter we confine ourselves primarily to the effect of time 
delayed coupling on the phenomenon of amplitude death 
and present a detailed numerical and analytical estimate 
of the parameter space in coupling strength and time de- 
lay where such a death can occur for identical oscillators. 
We also establish that this effect is not an artefact of the 
simple two oscillator model, but can occur for a system 
of large number of globally (or locally) coupled identical 
oscillators (including the continuum limit of iV ^ oo). 
We analyse the following model equations: 



Zi(<) = {i+iuji- 1 Zi(i) nzi(t) 

+K[Z2{t - t) - Zi{t)], 

Z2{t) - (l+^CJ2- I Z2{t) \^)Z2{t) 

+K[Zi{t - t) - Z2{t)], 



(1) 
(2) 



where r is a measure of the time delay, K is the cou- 
pling strength, aji_2 are the intrinsic frequencies of the 
two oscillators and Zi^2 are complex. The model is a 
generalization of the diffusively and linearly coupled os- 
cillators studied extensively for example in [p|pd|]. The 
time delay parameter is introduced in the argument of the 
couphng oscillator (e.g. Z2 in (1)) to physically account 
for the fact that its phase and amplitude information is 
received by oscillator Zi only after a finite time r (due to 
finite propagation speed effects). In the absence of cou- 
pling {K — 0) each oscillator has a stable limit cycle at 
I 1= 1 on which it moves at its natural frequency uoi. 
The coupled equations represent the interaction between 
two weakly nonlinear oscillators (that are near a Hopf bi- 
furcation) and whose coupling strength is comparable to 
the attraction of the limit cycles. It is important then to 
retain both the phase and amplitude response of the os- 
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dilators [||. The state Z^ = IS an eq 
of the system of equations [(Q) and (|2 



1 



equilibrium state is linearly unstable since the individual 
oscillators tend to stable limit cycle states | Zi |= 1. The 
stability of amplitude death for K ^ has been studied 
in great detail by Aronson et al [|j for system ((|l|,(|l)) in 
the absence of any time delay in the coupling (i.e. for 
T = 0). The conditions for stability found by them are, 

K>1 and A =| tJi - cjs |> 2^2^^, (3) 

which shows that amplitude death can occur in this case 
only for sufficiently large values of A provided K > \. 

In Fig. 1(a), we reproduce the bifurcation diagram 
of Aronson et.al ^ where the region marked / repre- 
sents the amplitude death region and the dotted curves 
mark the boundary as defined by condition (|^) . The two 
bounding curves intersect at the point [K — 1, A — 2). 
Regions marked {II) and (///) represent phase locked 
and phase drift regions respectively which we will not 
discuss in detail here. In Fig. 1(b) we present the bi- 
furcation diagram of (0-||) for r — 0.0817. Note that 
in contrast to the diagram of Fig. 1(a), the amplitude 
death region now extends down to A = and has a fi- 
nite extent along the coupling strength {K) axis. We 
find that the phenomenon persists for a range of t af- 
ter which the bifurcation curve lifts up from the A = 
line and identical oscillators can no longer suffer death. 
Fig. 2(a) shows this region, for different values of w, 
'm T — K space for which amplitude death of identical 
oscillators can occur. The size of this death island is a 
function of the frequency of the oscillators (w), as shown 
by the other curves. We shall soon show that the size 
is also a function of N , the number of oscillators. The 
bifurcation curves (including the island boundaries) have 
been obtained from a linear stability analysis of (0),(|^) 
about the origin {Zi = Z2 ~ Q) as well as direct numer- 
ical integration of the equations. Assuming the linear 
perturbations to vary as e^* the 
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FIG. 1. Bifurcation diagram of Eqns.(|l]-|). (a) t — 0. Re- 
gion / is the amplitude death region, // is the phase locked 
region, and /// corresponds to the phase drift or incoherent 
region, (b) r = 0.0817, and u) — 10. The death region ex- 
tends down to A = indicating that identical oscillators can 
suffer amplitude death. The phase locked region is split into 
two disjoint regions. 



characteristic eigenvalue equation we get is, 
{1-K + iuJi- A)(l - K + iuj2- X)- 

where A is the complex eigenvalue and the complete set 
of eigenvalues includes those arising from the complex 
conjugate equations of (1) and (2) . Setting Real{X) = 
in (y) and separating the real and imaginary parts, the 
equations for the critical curves (i.e. the marginal stabil- 
ity condition) are, 

Xj + 2XiQ + - — -{1-Kf 

+if^cos(2A/r) = 0, (5) 
2(1 - K){lu + Xi) + sm{2XiT) = 0. (6) 

where A/ = Imag(X) and u = (wi + W2)/2 is the mean 
frequency. Eliminating A/ between (|^) and (^) and con- 
sidering the full set of eigenvalues , we obtain the follow- 
ing transcedental relation between if, A and r which is 
now the modified marginal stability condition in place of 

ga — K'^ shi{aT ±2u)t)^ (7) 

where a = yA2~4g2"^4y^^^g^^^ j^^^j g = 1- K. 

Note that for r = 0, the above relation readily simplifies 
to (||) and yields the marginal stability curves K = 1 and 
2K = 1 -I- A_. Figure 1(b). is a numerical plot of (0) for 
r = 0.0817. 

To obtain a condition for the death of identical oscilla- 
tors we repeat the analysis with = a;2 = w in (^ and 
after eliminating A/, obtain the relations 
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FIG. 2. (a) The region of amplitude death for TV = 2 as a 
function of the common intrinsic frequency. The size of the 
the island decreases with decreasing frequency and vanishes 
below a certain threshold, (b) The death islands as a function 
of the number of globally coupled oscillators. Each oscillator 
is assumed to have an intrinsic frequency o{ u) — 10. The 
death island survives even in the limit of TV ^ 00. 
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The region of intersection between the two curves (||) 
corresponds to the death island region of Fig. 2(a). It 
demonstrates the only stabihty switches (of the origin) 
that take place as a function of r. For a given value 
of ll! and at a fixed K we move (by varying t) from an 
unstable region into a stable region as we cross the left 
boundary of the island to emerge again into an unstable 
region as we cross the right boundary of the island. No 
other amplitude death islands are seen for larger values 
of T (for a fixed w). We have also confirmed this ana- 
lytically by looking at the behaviour of the suprcmum of 
the real parts of the roots of the transcedental equation, 
dXii/dr (obtained from (4) with A = 0) as a function 
of r in the various parameter regimes . The detailed 
mathematical proof of this result will be published else- 
where. 

Can this phenomenon occur for an arbitrary number 
of oscillators? To answer this question we have inves- 
tigated the following generalized set of globally coupled 
equations: 



Z,{t) = [l + ioj,- I z,{t) V)z,{t) + 

N 

.Y^[Z,{t-T)-Z,{t)]- 
2K, 



2K 



N 



[Z^{t- 



(9) 



where i = 1,....,N and the last term on the right hand 
side has been included to remove the self-coupling term. 
For r = 0, (^) reduces to the set of equations that have 
been extensively studied by Ermentrout |9| , MiroUo and 
Strogatz §], and others |Q. MiroUo and Strogatz §] 
have provided rigorous analytical and numerical condi- 
tions for amplitude death in such a system. Their con- 
clusions, in general, are similar to the case of = 2, 
namely, that one needs a sufficiently large variance in fre- 
quencies for death to occur and K has to be sufficiently 
large. We have been able to carry out a similar linear 
stability analysis of (^ around the origin for the case 
of finite T and for a large number of identical oscillators 
{ujj — to, j — 1, N). The resulting stability condition 
yields the following bounding curves for the death island 



region: 



cos p 



oj - ^/AKh~l 
cos-H^p) ^ 



2tt — cos ^ p 
2n-cos-\j^p) 



(10) 



where p = 1 - l/{2Kb), q ^ AK^ ~ 1 + 4i^&(l - 2K) and 
the factor 6 = (1 — 1/A^) introduces the N dependence 
of the island size explicitly. In Fig. 2(b) we have plotted 
these islands for N = 2,3,4 and N — oo. To confirm 
these results we have also numerically scanned the re- 
gion with a direct numerical integration of (P) for a large 



number of oscillators, upto N = 10000, and found excel- 
lent agreement. We have also carried out a similar study 
for a large number of locally coupled identical oscillators 
(nearest neighbour coupling [||) with periodic boundaries 
and find that time delay introduces death islands in such 
systems as well. Thus it appears that in the presence of 
finite time delay in the mutual coupling, amplitude death 
of identical oscillators is a fairly universal phenomenon 
and occurs for any arbitrary number N of oscillators ex- 
tending upto N = oo over a range of t and K values. 
To the best of our knowledge such a result has not been 
realised in the past and may have important applications 
in biological or physical systems. There are many phys- 
ical examples of amplitude death in real systems. One 
of the earliest that was investigated both theoretically 
and experimentally is that of coupled chemical oscilla- 
tor systems e.g. coupled Belousov-Zhabotinskii reactions 
carried out in coupled stirred tank reactors |l^,|l^ . They 
can also occur in ecological contexts where one can imag- 
ine two sites each having the same predator-prey mech- 
anism which causes the number density of the species to 
oscillate. If the species are capable of moving from site 
to site at a proper rate (appropriate coupling strength) 
the two sites may become stable (stop oscillating) and 
acquire constant populations. Another important appli- 
cation of this concept is in pathologies of biological os- 
cillator networks e.g. an assembly of cardiac pacemaker 
cells jl^. Amplitude death signifies cessation of rhyth- 
micity in such a system which is otherwise normally spon- 
taneously rhythmic for other choices of parameters. For 
the onset of such an arrhythmia, current models based 
on coupled oscillator networks need to assume a signifi- 
cant spread in the natural frequencies of the constituent 
cells (oscillators) Our work demonstrates that this 
assumption may not be necessary if one takes into ac- 
count time delay effects arising naturally from the finite 
propagation times of the signals exchanged between the 
cells. Another possible application is in the area of high 
power microwave sources where it is proposed to enhance 
the microwave power production by phase locking a large 
number of sources such as relativistic magnetrons p2| . 
Time delay effects, arising from the finite propagation 
time of information signals traveling through the con- 
necting waveguide bridges, could impose important limi- 
tations on the connector lengths and geometries in these 
schemes. Our findings could provide a guideline in this 
direction. It should be noted that a form of oscillator 
death described in for identical oscillators is not a 
genuine amplitude death since it occurs in the context of 
a phase only model. Time delay in our study provides a 
new mechanism for genuine amplitude death to occur in 
coupled identical oscillators. 

Finally, it is worthwhile to mention that time delay can 
introduce other interesting phenomena as well, some of 
which have been studied in the context of the phase only 
model and need to be investigated for the more general 
phase and amplitude model. Our numerical results, for 
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FIG. 3. The bifurcation diagram of Eqs. (0) and for 
T — 0.4084 and Cj = 10. The amphtude death region (I) is 
surrounded by the phase locked regions (II) . The dotted curve 
which separates the incoherent (III) and the phase locked re- 
gions is obtained from numerical integration of the original 
equations. 

example, show that the bifurcation diagram of the sys- 
tem in the presence of time delay has a significantly richer 
structure. Fig. 3 is an example for the N = 2 system for 
T = 0.4084, which can be contrasted with the Aronson 
et al 1^ diagram of Fig. 1(a). Note that one no longer 
has the clean separation of the Bar-Eli region, the phase 
locked region and the phase drift region into three dis- 
joint regions that converge at a single degenerate point. 
Instead the phase locked region now always surrounds 
the Bar-Eli region and the single degenerate point is re- 
placed by a series of X points resulting from the braided 
structure of the phase locked region in the vertical direc- 
tion. At large values of K other bifurcation curves appear 
in the phase locked region indicating the appearance of 
higher frequency states A detailed investigation of 
various properties of this rich phase diagram, including 
stability studies of the various states, is now in progress 
and will be reported elsewhere. 
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